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We have used the XAFS (x-ray-absorption fine structure) technique to investigate the local 
structure about the Cu, Ba, and Tl atoms in orthorhombic Tl2Ba2CuOj, with a superconducting 
transition temperature T c =60 K. Our results clearly show that the O(l), 0(2), Cu, and Ba atoms 
are at their ideal sites as given by the diffraction measurements, while the Tl and 0(3) atoms 
are more disordered than suggested by the average crystal structure. The Tl-Tl distance at 3.5 
A between the TIO layers does not change, but the Tl-Tl distance at 3.9 A within the TIO layer 
is not observed and the Tl-Ba and Ba-Tl peaks are very broad. The shorter Tl-0(3) distance in 
the TIO layer is about 2.33 A, significantly shorter than the distance calculated with both the Tl 
and 0(3) atoms at their ideal 4e sites ( x = y =0 or i). A model based on these results shows 
that the Tl atom is displaced along the < 110 > directions from its ideal site by about 0.11 A; the 
displacements of neighboring Tl atoms are correlated. The 0(3) atom is shifted from the 4e site 
by about 0.53 A roughly along the < 100 > directions. A comparison of the Tl Lj/i-edge XAFS 
spectra from three samples, with T c =60 K, 76 K, and 89 K, shows that the O environment around 
the Tl atom is sensitive to T c while the Tl local displacement is insensitive to T c and the structural 
symmetry. These conclusions are compared with other experimental results and the implications 
for charge transfer and superconductivity are discussed. 

PACS numbers: 74.70.Vy, 74.75.+t, 78.70.Dm, 61.70.-r 



I. INTRODUCTION 



The Tl m Ba 2 Ca n _iCu n 2 (ri+i)+ TO system is a model family of high-temperature superconductors. The structures 
consist of copper perovskite-like blocks containing n Cu0 2 planes (n up to 6 has been observed) separated by one 
or two TIO layers. Each CaO layer is sandwiched between two Cu0 2 layers while the BaO layer is sandwiched 
between the Cu0 2 and TIO layers. For Tl-monolayer (m =1) phases, the superconducting transition temperature 
(T c ) increases to 122 K with the number of Cu0 2 planes from n — 1 to n — 4 and then decreases with n —5 to 6; 
for Tl-bilayer phases, T c increases up to 125 K with the number of Cu0 2 planes n from 1 to 3 and then decreases 
with n=4, 5. The fully stoichiometric Tl-bilayer Tl 2 Ba 2 Ca„_iCu ra 2rl +4 compounds would contain no carriers if 
the formal ionic valence states, Tl^+Baa+Ca^Cu^+O^, are assigned, and thus would not be superconductors. 
However, there is a hybridization of Tl 6s and O 2p orbitals which either touch the Fermi surface or cross it jj],^]], 
so the Tl does not act as a completely 3+ ion and charge transfer might occur between the Cu0 2 and TIO layers. 
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Other effects can also change the hole-doping in the system including: (i) excess or missing O; (ii) Ca 2+ and/or Cu 2+ 
substitution for Tl 3+ ; (hi) Tl deficiency; (iv) ordering of the O and/or Tl atoms in the TIO layers. Which doping 
mechanism primarily controls T c is still an unsolved problem. Structural studies on these compounds will be helpful 
in understanding the mechanism of superconductivity in this system. 

Tl2Ba2CuO y (Tl-2201) || is particularly suitable for this study because it has a simple structure (calcium-free and 
only 1 Cu02 layer per formula); the highest T c is about 90 K. The structure of Tl-2201, in fact, has been extensively 
studied since its first discovery in 1988 ||. Reasonable agreement for the general structure has been achieved from 
diffraction studies p|-|l2||. It contains TI2O2 double layers and a Cu02 plane separated by a BaO layer, as shown in Fig. 
1. For convenience of description, we define the 4e site {x = y = or ^) as the ideal site for Tl and 0(3) in tetragonal 
Tl2Ba2CuO y and treat orthorhombic Tl2Ba2CuOj, to be a pseudo-tetragonal structure. The atomic positions of Tl 
and O in the TI2O2 double layers, however, are not well defined, due to a large anisotropic thermal factor and/or 
various static displacements from their ideal sites. Most of the refined structures suggest a distorted site for the 0(3) 
with a displacement of ~0.31-0.42 A along the < 100 > direction or -0.22-0.28 A along the < 110 > 

direction ^,0,^), while some of them show that the Tl atom is also displaced from its ideal site by ~0. 11-0. 14 A along 
the < 110 > direction |7]j9| or —0.16 A along the < 100 > direction JTo| . The refined thermal factors in these studies 
are usually quite large in the ab-plane for the Tl and 0(3) atoms. Recently, it has been made clear that the thallium 
content has a strong influence on the lattice symmetry and thus on the average atomic position in the TIO layer 
jl2] Q Furthermore, significant distortions have recently been observed in both Tl2Ba2CuOj, and Tl2Ba2CaCu20s 
(Tl-2212) by means of x-ray-absorption fine structure (XAFS) |l5f| . Similar distortions are found in Tl-2212 by the 
pair-distribution-function (PDF) analysis of neutron scattering data [ ^6|Jl7| . It appears that these distortions not 
only have short-range correlations, but may in some cases also have long-range correlations, as evidenced by the 
presence of the modulational structure for selected sections, observed in high-resolution electron microscopy and 
electron diffraction |l8]-|2"2j] . In addition, the temperature dependence of the local structure in the Cu02 plane has 
been studied by XAFS in Tl-2201 jHll and h V PDF in Tl-2212 ||, respectively. The splitting of the Cu-O(axial) 
in TlBa2Ca3Cu40n has also been investigated by XAFS in the range 10<T<156 K p6[ . 

In this study, we have focused on the static distortion of the Tl2Ba2CuO l/ , especially in the TIO layer. We have used 
the XAFS technique |2^j2^| to probe the local structure around all the metal atoms, Cu, Ba, and Tl in Tl2Ba2CuOj,, 
to revise the earlier results based on the Tl L/77 edge studies []l5| . Both the sample preparation and the XAFS data 
analysis are slightly different from the previous study |ll|; they are described in Sec. II. The XAFS data analysis 
(presented in Sec. Ill) uses constrained fits |2t| to reduce the number of independent parameters and to check the 
previous model. The XAFS data and the data analysis for the Cu and Ba K-edges, and the Tl L///-edge are discussed 
in Sees. IV and V, respectively. A discussion and a comparison of our results with other experimental observations 
are given in Sec. VI, and the results summarized in Sec. VII. 



II. EXPERIMENTAL DETAILS 
A. Sample Preparation and Characterization 

Samples of about 10 grams weight were prepared from a mixture of TI2O3, Ba(OH)2-H20 and Cu(N03)2-nH20, 
always in the ratio of the cations desired in the final material with no excess of any element. The powders, milled 
in cthanol and thereafter dried in an inert atmosphere (N2), were heat treated by a two-step sintering process. The 
heating rate in both steps was 180° C/hour and the temperature was regulated by a controller, Eurotherm 818P, with 
a temperature stability better than ±0.5° C. 

In sintering step I the mixtures were heated in AI2O3 crucibles with AI2O3 covers at 760° C for 3 hours with only 
a limited loss, if any at all, of thallium. 

In sintering step II pellets were placed in AI2O3 crucibles sealed with gold gaskets and cover. They were subsequently 
centered in small pieces of ceramic tube and fixed with steel wedges (a procedure that prevents evaporation from the 
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crucible) and heated at 860° C for 6 hours. The loss of weight in sintering step II, if attributed solely to thallium 
evaporation, is below 0.3% of the total amount of thallium. There is no sign of substitution of gold into the material 
or of incorporation of thallium into the gold gasket. From EDX (energy dispersive x-ray spectroscopy) analysis, no 
traces of Al or Au can be detected in the sample. For estimated errors in barium and copper content see Ref. 12. 

Sample purity and lattice parameters were determined from Guinier film data using silicon (NBS 640b) as an 
internal standard. A computerized photoscanning system in combination with programs for indexing and unit cell 
refinements were used for the extraction of lattice parameters (a=5.4881(2) A, b=5.4532 A, c=23. 216(1) A). The 
impurity content was below 2%. 

Neutron powder diffraction data were collected at the pulsed neutron source ISIS, RAL, UK using the high flux, 
medium resolution powder diffractometer POLARIS. A least-squares profile refinement program, TF12LS, based on 
the Rietveld technique was used to analyze the TOF (time-of-flight) neutron powder diffraction data. The refinement 
was carried out in space group Abma and selected bond distances are found in Tables I-III. 

The superconducting properties were studied by ac susceptibility measurements using a Quantum design SQUID 
magnetometer at temperatures down to 4 K in a field of 50 Gauss. T c was determined to be 60 K. For further details 
concerning sample preparation and characterization consult Ref. 12. 



B. X-ray absorption measurements 

The XAFS samples were prepared by brushing a fine powder (w 30 /im) onto Scotch tape. Each layer was selected 
to be free of pinholes to minimize the XAFS glitches |30| and several layers were stacked to obtain a sample with 
a thickness of approximately one absorption length. Samples were measured at 80 K using a liquid-nitrogen dewar. 
Several data sets were collected for each sample to check for consistency and averaged to improve the signal-to-noise 
ratio. The x-ray absorption data were collected on beamline 10-2 at the Stanford Synchrotron Radiation Laboratory 
using Si (400) monochromator crystals. The double monochromator was detuned on its rocking curve to approximately 
one half of maximum transmitted x-ray intensity to reduce the harmonic content of the beam. 



III. XAFS DATA ANALYSIS TECHNIQUE 
A. Standard Procedures 

The standard procedures for the reduction and analysis of the XAFS data are described elsewhere ]27],^9|]. The 
extracted XAFS function, k\(k), is given by 

N F (k) 

k x (k) = 22 J ^ J 2 sin[2kRj + cj>j(k)]exp(-2k 2 a^ - 2R j /X) (1) 
j i 

where the sum is taken over shells with Nj atoms at a distance Rj from the absorbing atom, k is the momentum of 
the photoelectron given by k = [2m(E - Eo)} 1 / 2 /% where E a is the energy at one-half the absorption edge height, 
Fj(k) is the back-scattering amplitude, <j>j(k) is a phase shift of the photoelectron due to its interaction with the 
back scattering and absorbing atoms, A is the effective electron mean free path, and tjj (Debye- Waller factor) is the 
average variation (both static and dynamic) of Rj. Examples of the k-space data are given in Fig. 2 for the Cu and 
Ba K-edges and Tl L///-edge. A Fourier transform of these k-space data exhibits peaks in r-space which correspond 
to neighbors at different radial distances from the absorbing atom. 

To obtain numerical values for Rj , Nj and <jj , we perform iterated least-square fits of the real and imaginary parts 
of the Fourier transform of kx{k) to a sum of standards over a particular range in real (r) space. Identical k-space 
transforms are used for the sample and the standards. In these fits, the overlaps of the radial distributions of other 
neighbors in r-space are included. 
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It is difficult to extract all the isolated, single pair standards needed for this analysis from other experimental 
standards. In this study all the standards are generated by the FEFF5 code |3l]], using a cluster as close to the real 
lattice as possible with a radius of about 7.5 A from the central absorbing atom to the boundary scattering atoms. 
The multiple scattering paths are included in the fit when it is necessary. In earlier studies we have used experimental 
atom-pair standards extracted from data collected for standard compounds. We have compared our experimental 
standards with FEFF5 and have obtained very good agreement for about twenty pairs |32|. Since a complete set of 
experimental standards was not available, we decided, based on our recent study, to use FEFF5 standards. One of 
the problems in using the standards from FEFF5 is the uncertainty of the amplitude-reduction- factor, S%. We use 
two different values of S% which we have obtained from our detailed comparisons |52|; 0.8 when the back-scattering 
atom is O, and 0.9 when the back-scattering atom is a metal atom. 



B. Background- feature removal 

The standard background removal procedure which uses a simple spline or polynomial for fi (E), does not work well 
when the background consists of some features due to multielectron excitations |53],|4| and Ramsauer-Townsend-likc 
effects HQ. These features have been found at the Br (in RbBr), Kr, Rb (in RbBr), Ba (in BaO), Ce (in Ce0 2 ), 
and Pr (in PrBa 2 Cu30 7 ) K-edges §§-§§ and the Au, Tl (in T1 2 3 ), Pb (in /3-Pb0 2 ), and Bi (BaBi0 3 ) L/^-edges 
p2| , |33| in our studies. Similar features have also been observed and confirmed by other groups |3E| |. An iterative 
background extraction and removal procedure has been proposed to obtain and remove these features in the previous 
work §|. 

When the simple standard background removal procedure was applied to the Ba K-edge and Tl L/77-edge XAFS 
data of Tl 2 Ba 2 CuO a , spurious humps appear in the low-r region of the Fourier transformed data as shown in Fig. 
3. Instead of performing the iterative background removal procedure to these complicated data directly, we utilized 



the extracted background from the BaO |35(] and Tl 2 03 |32| data; the new background was obtained by fitting the 



corresponding extracted background plus simple spline functions to the data. These new backgrounds are plotted in 
Fig. 3. They contain more wiggles than the background obtain from the standard procedure but result in smoother 
r-space data in the low-r region than those from the standard procedure. This simple method works very well. We 
note that using more spline sections in the standard procedure may also minimize the humps in the low-r region but 
that usually reduces the amplitude of the first r-space peak, simultaneously. 



C. Constrained Fits 



One of the common problems in the XAFS data analysis is that there are too many fitting parameters if the fit 
includes further shells in a complicated system. For example, if the fitting range is 3.0 - 13.0 A" 1 in k-space and 1.0 
-4.0 A in r-space, the estimated allowed number of fitting parameters is ~21, according to the formula 2AfcAr/7r + 2 
p9[ . In a complicated system, such as the Cu edge data in YBa 2 Cu30y [ p9pc| , |4l| ] , the number of atomic pairs 
could exceed 8, which needs 24 or more fitting parameters, even if all the E c ' are fixed and not considered as fitting 
parameters. Therefore, a constrained fit is required in the XAFS data analysis for a complicated system. These 
constraints include using the lattice constants and the known average locations of the atoms in the unit cell from 
diffraction measurements, combining the XAFS results viewed by different atoms, and using equations which relate a 
parameter to one or more other parameters (for examples, keeping the ratio of some parameters constant; setting a 
parameter to the weighted sum of other parameters) to check possible models. The constraints force the fit towards 
a particular minimum in parameter space, and the fit is not easily trapped in some other local minima. In this way, 
the constrained fits can help reject results which are physically inconsistent and to check particular models for the 
local structure. This method has been used to a limited extent in our earlier work [fh5| , [4ll and more extensively in our 
recent studies [p9| . 
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IV. XAFS DATA PRESENTATION AND QUALITATIVE FEATURES 



A. Cu environment 

From diffraction studies [p-12|, the average local environment about a Cu atom is as follows: 4 0(1) neighbors 
at r=1.93 A in the CuC>2 plane, 2 0(2) neighbors at 2.71 A along the c-axis, 8 Ba neighbors at 3.35 A above and 
below the Cu02 plane, and 4 Cu neighbors at 3.87 A in the Cu02 plane. In XAFS data, the Fourier transform of the 
XAFS function k\{k) to real (r) space is a sum of peaks corresponding to each of these different neighbors; however, 
each distance is the average local value, not necessarily the value calculated from the average unit cell. The question 
considered in this section is how well does the Cu local environment agree with diffraction results. 

The Cu K-edge XAFS data in r-space are shown in Fig. 4. The transform range is a square window from 3.1 
to 16.5 A -1 , Gaussian broadened by 0.3 A -1 . Both the real part and the magnitude of the complex transform arc 
plotted. The distances between the absorbing atom and the back-scattering atoms are determined by the phase of 
the XAFS function, which is well represented by the positions of the zero-crossing points of the real and imaginary 
parts of the Fourier transform. However, there is a shift, —0.2 <~ —0.5 A, of the XAFS peaks from the real positions 
due the photoelectron phase-shifts. 

To compare the local structure which determines the Cu K-edge XAFS with the average structure obtained from 
diffraction, we generated the theoretical XAFS function using FEFF5 |5l| and the refined crystal structure obtained 
from diffraction measurements Q|. In the calculation, all the atoms within a distance of 7.5 A from the central 
absorbing atom, and all the scattering paths with a contribution of 4% or more, relative to the maximum contribution, 
are included. To roughly simulate both the static and thermal disorder in the calculation, we have set the sample 
temperature to 100 K and the Debye temperature to 700(K A)/(0.5Rj) which increases a with distance Ft, for more 
distant shells. (The same thermal parameters are used for the other edges in this paper unless noted). The Fourier 
transform of the theoretical XAFS (without any fitting) is shown in Fig. 4 by dotted lines, with the same k-space 
window as for the experimental data, an overall amplitude-reduction-factor of 0.9, and an overall E Q shift of —2.0 eV. 

Fig. 4 shows that the phase of the theoretical XAFS agrees very well with that of the experimental data up to 
6.5 A. This means that the 0(1), 0(2), Ba, and Cu atoms are well ordered about a given Cu atom, in excellent 
agreement with the diffraction results. We point out that the Cu edge XAFS data is not sensitive to the distortion 
of the distant TIO layers if there is any, especially when the distortion is in the ab-plane, because the TIO and Cu02 
layers are separated by 4.7 A and a small splitting of the Cu-Tl and/or 0(3) distances can be absorbed by a large 
Debye- Waller factor. We note that the subtle differences in zero-crossing points are due mainly to the errors in E D , 
especially in the region below 2.3 A where the O atoms make a major contribution. The nearest O neighbors directly 
touch and partially screen the excited absorbing atom which has a hole in the core-shell level, thus the E Q shift for the 
nearest O neighbors should be different from the distant neighbors. Since the same E Q shift is used for all neighbors 
in the simulation, an error in the phase-shift can be introduced in the first O peak. The discrepancy in the magnitude 
is mainly from the rough estimation of the Debye- Waller factors and the amplitude-reduction-factors of the XAFS 
function in the theoretical calculation and is not strongly related to the separations of the atomic pairs. 

We emphasize that in the theoretical calculation we did not adjust the XAFS parameters, Nj, Rj, Uj, and E OJ -, 
individually; instead, we only adjusted the overall amplitude and E G and simulated the Debye temperature of the 
different pairs by a simple formula. This method makes it possible and easy to compare the experimental XAFS data 
with a theoretical function calculated for a particular model up to 7 A (for example, from the refined average unit cell 
obtained from diffraction). However, it is still difficult to make a detailed fit to the experimental XAFS data much 
beyond 4 A. 



B. Ba environment 



The Ba atom is in a plane between the Cu02 and TIO layers. From diffraction measurements NH12J], it has 4 0(1) 
neighbors at 2.73 A, 4 0(2) neig hbors at 2.84 A, 1 0(3) neighbor at 3.00 A, 4 Cu neighbors at 3.35 A, 5 Ba neighbors 
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(4 in the ab-plane and 1 along the c-axis) near 3.86 A, and 4 Tl around 3.9 A 0). Again we ask the question if the 
local structure about the Ba atoms is the same as that given by the diffraction results assuming no distortions of the 
Tl atoms. The Fourier transform of k\(k) for the Ba K-edge data is shown in Fig. 5; window is from 3.5 A^ 1 to 15.5 
A -1 , Gaussian broadened by 0.3 A -1 . We have generated the theoretical Ba XAFS function following the procedure 
described earlier for the Cu edge; it is plotted as dotted lines in Fig. 5a. 

The phase of the theoretical and experimental XAFS function agreed very well up to 3.7 A (a small discrepancy in 
the phase of the first peak can be easily minimized by using a different E D shift). This again indicates that the 0(1), 
0(2), Cu and Ba atoms are well ordered. Similarly, there is excellent agreement in the phase for a large peak around 
5.2 A between the simulation and experimental XAFS data. The dominant contributions to this peak arise from 
the collinear, forward, multiple scattering of the Ba-0(1)-Ba and the multiple scattering of intraplanar Ba-0(2)-Ba 
pairs. The phase of the real part of the transform within the 5.2 A peak is shifted by nearly 7r/2 relative to the phase 
that would be observed for a single scattering Ba-Ba peak (e.g. as occurs near 3.5 A in Fig. 5) due to the forward 
multiple scattering contributions. The good agreement between the calculations and the experimental data for this 
peak strongly suggests that the Ba atom is well ordered. 

More importantly, however, there is a significant difference between the calculations and the experimental data on 
the high-r side of the peak near 3.6 A(see Fig. 5c in which a small fraction of the data is shown with an enlarged 
x— axis), where the Ba-Tl pair (r= 3.90 A) should make a major contribution. The experimental data in this region 
of r-space can be simulated much better by FEFF5 when the Ba-Tl pairs are excluded; i.e., the net Ba-Tl peak is 
small. This means that the Ba-Tl peak is either very broad due to disorder of the Tl atoms or is split into multi-peaks 
due to a displacement of the Tl atom from its ideal site. 

C. Tl environment 

There are two adjacent T10 layers, which are surrounded by two BaO layers, one above and one below. The average 
local structure around the Tl atoms obtained from diffraction [Q-[l2| consists of 1 0(2) and 1 0(3) atoms near 2.0 A, 
4 0(3) within 2.3-3.2 A; 4 interplanar Tl atoms near 3.5 A and 4 intraplanar Tl atoms near 3.9 A; and 4 Ba atoms 
near 3.9 A. However as outlined in the introduction and the above discussion at the Ba edge, there is considerable 
evidence that the T10 layers are distorted. Here we investigate this disorder using Tl edge XAFS data. 

First, we performed a multiple-scattering-cluster-simulation for the Tl L/77-edge XAFS (as described above for the 
Cu K-edge data) assuming the structure of Ref. 4 in which the Tl is at the ideal 4e site while the 0(3) is displaced 
along the < 100 > direction by 0.37 A from its ideal 4e site. Both the simulated and experimental XAFS, k\{k), were 
Fourier transformed to r-space using a k-space range of 3.2 - 15.5 A -1 , Gaussian broadened by 0.3 A -1 , as shown 
in Fig. 6a. In the simulation, we used two different sets of Debye- Waller temperatures, 300 (K A)/(0.5Rj) for Tl-Tl 
pairs and 700 (K A)/(0.5i?j) for other pairs, to better simulate both the static and thermal disorder. This gives a 
good simulation above 4.5 A in r-space. 

Significant differences in the phase between the simulated and experimental data are observed in Fig. 6a. Thus at 
least some of the distances derived from the average, crystallographic structure are not consistent with these XAFS 
data. For the first O neighbor peak (1.3-2.3 A in Fig. 6) of the simulated data, the phase of the low-r side (1.3-1.8 
A), where the Tl-0(2) and Tl-0(3) pairs at r=2.0 A make their major contributions, agrees reasonably well with the 
experimental data; however, at the high-r side of this peak (1.8-2.3 A) , where the Tl-0(3) pair at r=2.5 A makes 
its major contribution, the phase is not consistent with the experimental data and a contribution in the simulated 
data is missing in the experimental data. A comparable inconsistency is observed in the peak centered near 3.6 A. 
The phase and the shape of the simulation for this peak agree well with the actual data in the low-r side, where 
the interplanar Tl-Tl pair makes its major contributions; however, it does not agree in the high-r side (3.7-4.0 A), 
where the intraplanar Tl-Tl and also the Tl-Ba pairs make their major contributions. Furthermore, the amplitude 
of the simulated peak in r-space near 3.6 A is almost twice as big as that in the data. This discrepancy is too large 
to be explained just by the errors in the amplitude-reduction-factor and the Debye- Waller factors. This implies a 
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significant distortion and/or disorder in the Tl-Tl and Tl-Ba distances from the average distances obtained from the 
refined structure in diffraction. A simulation without the intraplanar Tl-Tl pairs bears a much closer resemblance to 
the experimental spectra near 3.5 A. These comparisons again suggest that the Tl atoms, and probably also the 0(3) 
atoms, occupy more distorted sites in the ab-plane than those obtained from the average diffraction experiments. 

V. DETAILED XAFS DATA ANALYSIS 

From the qualitative analyses of the Cu, Ba, and Tl edge XAFS data shown above, we conclude that the positions of 
the Cu, 0(1), 0(2), and Ba atoms in the unit cell are consistent with the average crystallographic structure obtained 
from diffraction studies. However, the Tl and 0(3) atoms are significantly more distorted than the average structure 
obtained in many diffraction studies. The XAFS data indicate that these distortions are located in the ab-plane of 
the T10 layer as suggested in some diffraction investigations. These qualitative features are model independent and 
give a sound guidance for further structural analysis. 

In this section, we will perform nonlinear-least-square fits to the XAFS data to further check the qualitative features 
shown above and gain quantitative information on the local structure about the Cu, Ba, and Tl atoms. 

A. Cu environment 

The r-space Cu K-edge XAFS data have been fit up to 3.9 A using individual standard atomic pairs generated by 
FEFF5. Five pairs are included in the fit: Cu-O(l) near 1.9 A, Cu-0(2) near 2.7 A, Cu-O(l) near 4.3 A, Cu-Ba near 
3.4 A, and Cu-Cu near 3.9 A in which the multiple-scattering-paths for Cu-0(1)-Cu are included. The experimental 
data and the fit are shown in Fig. 7, where the same Fourier transform range is used as in Fig. 4 and the fitting range 
is between 1.0-3.9 A. Overall, a high quality fit has been achieved. A double-well potential was used in Ref. 23 to 
improve the fit for the first Cu-O(l) peak; however, in this study, the refinement of the first peak is not important and 
we use only a single Cu-O(l) peak. The structural parameters extracted from the XAFS curve-fitting are tabulated in 
Table I, together with the results from diffraction refinement. Our XAFS results are in excellent agreement with the 
diffraction results, as expected from our earlier comparison of the simulated and experimental XAFS r-space spectra, 
in which there was an excellent phase agreement up to 6.7 A. These XAFS results are also very close to the previous 
results (42). 

B. Ba environment 

In Sec. IV, we showed that the Tl is disordered or distorted. So in the fits to the Ba K-edge data, we considered 
three situations for comparison: (1) no Ba-Tl peak, (2) a single Ba-Tl peak, and (3) a split Ba-Tl peak. For the first 
case, we fit the Ba K-edge XAFS data with 3 O peaks (at r=2.7 A, 2.8 A and 3.0 A, respectively), 1 Cu and 1 Ba 
peak. The Tl peak near r=3.9 A was not included. The Debye- Waller factor of the Ba-Cu peak was fixed to have the 
same value as that of the Cu-Ba peak from the Cu K-edge fit. A reasonably good fit is achieved and all the extracted 
parameters are reasonable, which means that the Ba-Tl pairs do not make a significant contribution to the Ba edge 
XAFS, in contradiction to what one would expect for the undistorted Tl layers. 

Secondly, we included the Ba-Tl peak in the fit and fixed the amplitude of that peak to 4, the number predicted 
by diffraction. The quality of the fit was improved by about 25% and the extracted distances and the amplitudes for 
the Ba-O, Ba-Cu and Ba-Ba peaks remained essentially unchanged. The Ba-Tl distance is about 3.90 A, very much 
the same as the average distance given by diffraction. The Debye- Waller factor of the Ba-Tl peak is about 0.094 A, 
much bigger than that of the Ba-Ba peak (0.066 A). 

Finally, we tested the model proposed earlier |ljj . Here we require consistency between the XAFS spectra collected 
from three different absorbing atoms, achieved by using constraints in the fitting procedure. Following our previous 
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discussion, we assumed that the Tl atoms are shifted along < 110 > so that the 4 Ba-Tl pairs split into three distances, 
with roughly 1 pair at r — Sr, 2 pairs at r and 1 pair at r + Sr. Therefore, we included 3 0,1 Cu, 1 Ba and three 
Ba-Tl peaks in the fit. To minimize the number of fitting parameters and to check the model, we fixed the average 
distance r and the amplitudes of these three peaks to the values predicted by the model and only allow Sr to vary 
for these Ba-Tl peaks. The Debye- Waller factors of the three Ba-Tl peaks and the E Q shifts for the same type of 
back-scattering atoms were also kept equal. In this way, the number of the variables was reduced to 18, well below 23, 
the estimated number of allowed variables obtained from (2 Ar Ak/n + 2) |39] . The quality of the fit and the extracted 
parameters for O, Cu and Ba peaks are essentially the same as the fit with a broader Ba-Tl peak but here i5r=0.10 
A. As a final step, we also allowed the average Ba-Tl distance and the amplitudes of the Ba-Ba and Ba-Tl peaks to 
vary. The final fitting parameters remained essentially the same as the ones with the strong constraints. 

The experimental spectra and the final fitting curves are plotted in Fig. 8a; the contributions from O, Cu and Ba 
atoms are shown in Fig. 8b, while the contributions from the split Ba-Tl pairs are shown in Fig. 8c. These results 
clearly show that the contributions of the Ba-Tl pairs are quite small compared to the Ba-Cu and Ba-Ba pairs. The 
detailed results are listed in Table II, together with the results from diffraction and the model for comparison. We 
would like to point out that the Ba-0(3) distance is not sensitive to the displacement of the 0(3) atom in the ab-plane. 
For Ba-Tl peaks, Sr is slightly coupled with the Debye- Waller factor a and the coordination number N. It varies 
from 0.08 A to 0.11 A depending on a and N (or S%). We also used different k weighting (i.e., k n x(k) with n=l, 2, 
3) and different k ranges for the Fourier transform in the fit. The derived parameters are essentially the same within 
the estimated uncertainty. In all cases, the splitting distance, Sr, for the Ba-Tl pair is significantly shorter than 0.2 
A, the distance predicted by the previous model. This means the displacement of the Tl atom in the previous model 
(0.28 A) is too large and has to be modified to fit the Ba edge XAFS data. 

We also simulated the Ba K-edge XAFS data with a distorted model for the TIO layer. Here we modified the 
previous model by assuming that the Tl atom is displaced 0.11 A from its ideal site along the < 110 > direction. The 
simulated data are plotted in Fig. 5b. The agreement with the data is much improved for the r-space range 3.5-4.5 
A compared to the earlier simulation (Fig. 5a), as illustrated in Fig. 5c. There is also improved agreement in phase 
near 6 A. 



C. Tl environment 

From the analysis of the Cu and Ba K-edge XAFS spectra, we conclude that the 0(1), 0(2), Cu, and Ba atoms are 
at their ideal sites in the unit cell just as diffraction shows, but the Tl atom is shifted from its ideal site. However, it 
is difficult, from the Cu and Ba edge XAFS spectra only, to distinguish between a randomly disordered Tl site (which 
gives a broad peak) and a distorted Tl site (which gives multiple peaks) and hence to derive the actual displacement 
of the Tl atom if it is at a distorted site. Further information from the Tl edge is desirable. 

The Tl L//7-edge XAFS data are more sensitive to the Tl distortion, but the data analysis is more complicated by 
the fact that there are two adjacent TIO layers. In the fits to the Tl edge data, we considered the Tl displacements 
along the < 110 > direction with two different magnitudes: (1) ~0.3 A, similar to the previous model |15|, and (2) 



~o.i A, as suggested by the Ba edge XAFS data. In comparisons with diffraction studies, we also tried other models, 
see Sec. VI.A. 

In order to enable a meaningful analysis of the Tl L///-edge data we treat the Tl atoms as if they all see the same 
local distortions. However, we know that in the related Bi-based cuprate superconductors, there exist strong periodic 
modulations |^| which result in superstructure with long range order, a buckling in the ab-plane, and periodically 
extra O atoms. Because diffraction results show little disorder for Tl along the c-axis, we are motivated to use a 
simple model for the local environment. As discussed later, it is still compatible with a superstructure unit cell. 

We first performed a rough fit to the data in the region of 1.3^3.8 A, with 2 O, 1 Tl and 1 Ba peaks, similar to the 
model used in the previous paper p5[ . The extracted structural parameters are the following: 1.8 O atoms at 2.04 A, 
2.9 O atoms at 2.35 A, 2.8 (which is too small compared to 6) Ej Tl atoms at 3.51 A, and 7.6 Ba atoms (but very 
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broad, with a = 0.13A) at 3.89 A. This shows that the interplanar Tl-Tl pair (at r=3.5 A) remains but the intraplanar 
Tl-Tl pairs (with an average distance of 3.9 A) do not make a large contribution to the Tl edge XAFS spectra as 
would occur for an undistorted site. This result suggests that the Tl atom is distorted (not randomly disordered) in 
a particular way so that the interplanar Tl-Tl peak remains but the intraplanar Tl-Tl peak seems to disappear. The 
fitting result we obtain here is quite close to the earlier result M| . However, the model based on that result requires 
a large shift for the Tl atom which does not fit the Ba edge data as discussed above, and the Tl-Tl amplitude is too 
low. Therefore, the previous model must be modified to fit both the Ba and Tl edge XAFS data. 

Secondly, we followed the previous model but with a smaller shift for the Tl atom. Assuming that the Tl atoms 
are shifted along the < 110 > direction as shown in Fig. 9, the interplanar Tl-Tl distance remains approximately 
unchanged while the intraplanar Tl-Tl distance splits into two distances. In addition, the Tl-Ba distance splits into 
three distances, with the average distance unchanged to first order. In order to check this model more carefully, we 
explicitly included all the split peaks and the long Tl-0 peak at 4.3 A in the fit using the constrained-parameter 
technique. The numbers of neighbors for all the Tl and Ba peaks are fixed to the values predicted by the model. 
The average Tl-Ba and intraplanar Tl-Tl distances are also fixed to the average distances given by the diffraction 
measurement, and only the split distances are allowed to vary. The Debye- Waller factors of the three Tl-Ba peaks 
are set to be equal to those of the Ba-Tl peaks at the Ba edge, to further reduce the number of independent variables 
and to check the consistency of the fitting results. The E D shifts for the same type of back-scattering atoms are kept 
equal except for the nearest O atoms. Thus the total number of variables is 16, far below the estimated number of 
allowed variables, 27 |p9|| . The fit was carried out in the region of 1.2 A< r <4.2 A, with the Fourier transform range 
of 3.3 A -1 ~ 15.5 A -1 , Gaussian broadened by 0.3 A -1 . A high quality fit was achieved. 

Finally, we allowed the average Tl-Ba and intraplanar Tl-Tl distances, the amplitudes, and the Debye- Waller for 
the three Tl-Ba pairs to vary. The fitting results remained essentially the same as those with the constraints. We also 
fit the data with different k weighting (k n x(k), n=l, 2, 3) and k range in the Fourier transform. Consistent results 
were achieved within the estimated uncertainty. 

The final fit is shown in Fig. 10a. The contributions from O, Tl and Ba atoms are plotted in Fig. 10, b-d, 
respectively. The contributions from the three individual Tl-Tl peaks are plotted below the sum of their contributions. 
The contribution of the Tl-Tl at r=4.02 A is small, and the Tl-Tl peak at r=3.72 A is out of phase with the Tl-Tl 
peak at r=3.50 A, resulting in partial destructive interference. Thus the sum of the Tl-Tl peaks has almost the same 
phase as that of the Tl-Tl peak at r=3.50 A but a smaller amplitude, which explains why the fit with only one Tl-Tl 
peak gives quite good agreement, but with a small amplitude. The three Tl-Ba peaks are also plotted below the sum 
of their peaks. The sum of the three split Tl-Ba peaks is a broad peak, with a similar phase to that of the Tl-Ba peak 
at the average diffraction distance, 3.90 A, which again is consistent with one broad Tl-Ba peak in the previous fit. 

The extracted structural parameters are listed in Table III, together with the corresponding parameters calculated 
from the model. We emphasize that the Tl-Ba distances obtained in this new fit of the Tl L///-edge data are in good 
agreement with the Ba-Tl distances extracted from the Ba edge XAFS data within the estimated uncertainty. In 
the model, we only shifted the Tl atom along the < 110 > direction to match the Tl-Tl, Tl-Ba and Ba-Tl distances 
extracted from our XAFS data. The best overall agreement is achieved when the Tl atom is displaced by about 0.11 
A. We have also shifted the 0(3) atom along the < 100 > direction from the 4e site by 0.53 A in the model to match 
the observed average interlayer Tl-0 distance at 2.33 A. The splitting of the two shorter interlayer Tl-0 distances 
(r=2.29 and 2.38) predicted by the model is consistent with the broad Tl-0 peak (<r=0.131 A) at r=2.33 A in the 
data. We point out that absolute position for the 0(3) is not well defined in the data because the Tl-0 peak at 
2.33 A is quite broad and the Tl atom itself is displaced. The long interlayer Tl-0(3) peak, predicted to be near 
r =3.12-3.24 A by the model, is not clearly observed in the XAFS data, probably because that peak is considerably 
broader than the shorter one, and/or smeared out due to multiple Tl-0 peaks in that range. Instead of the well 
defined pair distances predicted by the model, the actual PDF of the Tl-0 (3) pairs within the TIO layer may be 
spread over the entire range 2.3-3.2 A (refer to Table III) due to the variations of the Tl and 0(3) displacements in 
both magnitudes and directions. 

We again simulated the Tl edge XAFS, using the FEFF5 code and the distorted local structural model (see Fig. 
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9). The simulation is shown in Fig. 6b. The sample temperature and Debye- Waller temperature used here are the 
same as those used in the simulation shown in Fig. 6a, except the parameters for the first two Tl-0 pairs (r=2.00 A 
and 2.09 A ) in which the Debye- Waller factors were set to 0.03 A because the static disorder had already partially 
been taken into account by the two separated peaks. The new simulation in Fig. 6b shows a great improvement over 
the old one shown in Fig. 6a, and the phase of the new simulation agrees very well with the actual data over the 
whole region. 

VI. DISCUSSION 
A. Comparisons with diffraction results 

The extracted distances from our XAFS data analyses for the undistorted atomic pairs, such as Cu-O, Cu-Ba 
(Ba-Cu), Cu-Cu, Ba-0(1,2), Ba-Ba, are in excellent agreement with those obtained from diffraction studies, with a 
typical error of 0.01 A (see Table I and II). This shows the high quality of the XAFS data and the data analyses. The 
extracted distances for the distorted atomic pairs, such as Ba-0(3), Ba-Tl (Tl-Ba), Tl-O, are within the range of the 
distances given by diffraction investigations (see Table II and III). 

Diffraction studies Pf^-plPI show that the refined thermal factors for the Tl and 0(3) atoms along the c-axis are 
only one third or less of those in the ab-plane. Therefore, we did not attempt to vary the c-components of the Tl and 
0(3) positions in our model; they are fixed at the average positions given by diffraction studies. Furthermore, the 
mean-square-root of the Tl displacement is about 0.15 A in the ab-plane (y/Uu or \/U%i) and 0.09 A along the c-axis 
(VU33) in the refinement of diffraction studies P,[7[-p|,pT[ , which are quantitatively consistent with our XAFS results. 
Given the Tl displacement of 0.11 A in the ab-plane and the Tl Debye- Waller factor of 0.09 A (see Table III), the 
approximately estimated mean-square-root of the Tl displacement in the ab-plane is 0.14 A (VO.ll 2 + 0.09 2 A), very 
colse to that of the diffraction results. 

The model (Fig. 9) based on our XAFS data analyses suggests that the Tl displacements are correlated over a 
short range. The partial-pair-distribution-function (PPDF) for the Tl-Tl pairs extracted from our data shows a dip 
near the average distance as shown in Fig. 11a, which is significantly different from the Tl-Tl PDF for the Tl in the 
ideal sites as shown as the dotted curve in Fig. lib (with a Debye- Waller factor 0.07 A). The second peak in the 
ideal structure is split into two peaks in the extracted PDF, one near 3.7 A and another one near 4.0 A, thus a dip is 
introduced near the average distance 3.87 A. When we fit the data using the ideal structural model (4 Tl-Tl at 3.51 
A, 4 Tl-Tl at 3.87 A and 4 Tl-Ba at 3.89 A), the longer Tl-Tl pair is moved to 3.75 A (solid curve in Fig. lib) which 
physically does not make sense because the average Tl-Tl distance in the ab-plane should equal to the average lattice, 
3.87 A. 

Our model is similar to that suggested by the PDF analysis of the pulsed-neutron scattering spectrum in Tl-2212 
0,0, but more specific about the relative shift for the Tl atoms between the two consecutive TIO layers. Since 
the interlayer Tl-Tl distance remains unchanged, the displacement of the nearest Tl-Tl pairs between the TIO layers 
should be either in the same direction or perpendicular to each other, but not in opposite directions. The shifts of 
the Tl and 0(3) atoms and the pair distances obtained from XAFS, however, are different from those obtained by 
PDF even if the small changes in the lattice constants between Tl-2201 and Tl-2212 are taken into account. The 
PDF analysis for Tl-2212 suggests that the Tl atom shifts by 0.32 A along the < 110 > direction, which results in 
three split Tl-Ba distances with Sr ~0.23 A. This is not consistent with our Ba K-edge XAFS data (Sr ~0.1 A). We 
point out that XAFS probes the local environment around each selected element separately, and thus determines the 
average local structure around every different element. The combination of the information from different types of 
absorbing atoms gives a more constrained result for the local structure. On the other hand, PDF analysis uses all the 
atomic pairs in the material together, so that identifying individual partial pair-distribution-functions (PPDFs) is a 
much more involved process. Consequently, XAFS offers an advantage in studying the near neighbor local structure 
around an absorbing atom in complicated materials. 
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Some diffraction studies also revealed a Tl displacement from its ideal site [|^p|JlO|]. An x-ray structural study of 
tetragonal single crystals (T^ 85 Cuo.i5)Ba2Cu06 with T c =110 K suggested that the Tl is displaced along the < 110 > 
direction by 0.11 A ||. Another x-ray diffraction study on tetragonal single crystals Tli.gyE^Cui.iaOg+a with T c =12 
K showed a Tl displacement of 0.16 A along the < 100 > direction po| . A neutron powder diffraction investigation 
of orthorhombic T^I^CuOe+a with T c =90 K inferred a Tl displacement of 0.14 A along the 6-axis (which is 
along the < 110 > direction in a tetragonal unit cell). In all these diffraction studies, the Tl atoms are assumed to 
be statistically distributed in the displaced sites. This leads to an uncorrelated Tl displacement with the Tl-Tl PDF 
having its largest amplitude at the average distance (see dotted lines in Fig. lib and c), which is significantly different 
from the PPDF proposed by our XAFS study. To further check the models suggested by diffraction studies, we have 
tried to fit the XAFS data with the uncorrelated Tl displacement model for displacements along both the < 110 > 
and < 100 > directions. 

1. Tl displacement along the < 110 > direction 

In the < 110 > direction Tl displacement model suggested by one diffraction study ||], the four- fold displaced sites 
of the Tl are randomly occupied. Thus, to check that model in the XAFS data analysis, eight Tl-Tl pairs must be 
included in the fit with distances r% — 25r\, r% — 5r\, 7*1, r\+ Sr%, r\ + 2<5t*i, r2 — 1.9Sri, r2, r2 + 1.9Sr± and amplitudes 
0.25, 1, 1.5, 1, 0.25, 1, 2, 1, respectively. Here, r\ and are the average Tl-Tl distances between and within the 
TIO layers, respectively; the splitting of the Tl-Tl distances from their average distances are proportional to Sri to 
first order. Correspondingly, three Tl-Ba pairs must be included in the fit with distances r^ + Sr^, r^, r% — Sr% and 
amplitudes 1, 2, 1, respectively. Three Tl-0 pairs near 2.0 A, 2.3 A and 4.4 A were also included. 

The fit was started with the Tl displacement of 0.11 A from its ideal site (which corresponds to <5ri=0.085 A and 
<Sr 3 =0.08 A) and Debye- Waller factors of 0.07 A for the Tl-Tl and Tl-Ba pairs. The starting Tl-Tl PDF is shown 
in Fig. lib as dotted line. To reduce the number of fitting parameters and to check the proposed model, we fixed 
n, 7*2, and r 3 at the average distances (3.51 A, 3.87 A, and 3.89 A, respectively) given by diffraction study, and 
the amplitudes of the Tl and Ba peaks to the numbers predicted by the model. The E D shifts for the same type of 
back-scattering atoms are kept equal except for the nearest O atoms. With these constraints, the total number of 
variables is 25 for this 14-pair-fit. The fitting range is exactly the same as that used for the correlated model (Sec. 
V.C). 

The extracted PPDF is plotted in Fig. 11c as the solid line. The quality of this fit with 25 parameters is comparable 
to the one obtained from the constrained fit using the correlated model (Fig. 9) with only 16 variables. The distances 
changed a little from the starting values to Sr\ =0.095 A and dr^—0.07 A. Since the weighted number of the Tl-Tl 
pairs are fixed, their Debye- Waller factors are varied to adjust the height of the peaks. The extracted Debye- Waller 
factors are less than 0.04 A for the Tl-Tl peak at 3.51 A and larger than 0.2 A for the Tl-Tl peak at 3.87 A, which 
results in a sharp peak at 3.51 A and a dip near 3.87 A in the Tl-Tl PDF as shown in Fig. 11c by the solid line. For 
the Tl-Ba PDF, there is no difference between the < 110 > direction correlated and uncorrelated Tl displacements, 
and the extracted Tl-Ba PDFs are very similar. 

Although the quality of the fit for the uncorrelated model (25 variables in the fit) is comparable to that of the 
correlated model (16 variables in the fit), the extracted Tl-Tl PDF (solid line in Fig. 11c) from the constrained fit is 
significantly different from the starting Tl-Tl PDF ( dotted line in Fig. 11c) which corresponds to the uncorrelated 
model. Furthermore, the extracted Tl-Tl PDF appears to be very close, within the imposed constraints, to that 
obtained from the correlated displacement model (solid line in Fig. 11a); a dip shows up at the average distance 3.87 
A. For the uncorrelated model, a peak must show up at the average distance in the Tl-Tl PDF. Clearly, our XAFS 
data do not support the uncorrelated < 110 > Tl displacement model. 
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2. Tl displacement along the < 100 > direction 



For the < 100 > direction uncorrelated, four-fold Tl displacement model suggested by diffraction study |10fl , eight 
Tl-Tl pairs must be included in the XAFS fit with distances r\ — 5ri, r±, r\ + Sri, r 2 — \.llSri, r 2 — 0.885<5ri, r 2 , 
r 2 + 0.885<5ri, r 2 + 1.77Sn and amplitudes 1.0, 2.0, 1.0, 0.25, 1.0, 1.5, 1.0, 0.25, respectively. Two split Tl-Ba pairs 
are included in the fit with distances r^ — Sr^, r% + Sr%. Three Tl-0 pairs are also included. 

The fit was started with the Tl displacement obtained from diffraction, 0.16 A |[o| (which corresponds to Sri = 0.15 
A and <5r 3 =0.07 A). As was done for the fit to the < 110 > uncorrelated Tl displacement model, ri, r 2 , r%, and the 
amplitudes of the Tl and Ba peaks are fixed. The total number of variables is also 25. The quality of the constrained 
fit is about 30% worse than that of the fit for the < 110 > direction correlated displacement. The Tl-Tl peak at 3.87 
A is depressed by a large Debye- Waller factor, 0.2 A, in the fitting process. The obtained distances, with Sri— 0.25 A 
and 5r3=0.06 A, deviate significantly from the starting model and are no longer consistent with the starting model. 
The extracted Tl-Tl PDF is plotted in Fig. lid (solid line), which clearly moves away from the starting Tl-Tl PDF 
suggested by diffraction and towards the correlated Tl-Tl PDF (solid line in Fig. 11a). 

The extracted Tl-Ba PDF is somewhat similar to that obtained from the < 110 > direction displacement, but it is 
difficult to distinguish between one broad peak, two split peaks, and three split peaks. 

The above analyses clearly show that the Tl-Tl peak at the average lattice distance, 3.87 A, is essentially destroyed 
by the large Debye- Waller factor in the uncorrelated Tl displacement model, for displacements along either the < 110 > 
or the < 100 > directions, and the extracted PPDF from the uncorrelated model tends to move away from the PPDF 
generated by the starting model and towards the PPDF obtained by the correlated model. Therefore, we conclude that 
our XAFS data do not support the uncorrelated Tl displacement along either the < 110 > direction or the < 100 > 
direction proposed by the diffraction studies |7|,P,^0|. We point out that XAFS is more sensitive than diffraction 
to small local structural distortions which have no long range order. This can be seen from the following example. 
Rietveld refinements were performed for Tl-2212 using 7 distorted models with different directions of displacement 
including two correlated PDF models jl7|| . The results of these refinements show relatively little sensitivity to the 
choice of model used for the TIO layers. Three models, i.e., the uncorrelated four-fold displacement of both 0(3) and 
Tl along the < 100 > directions, the uncorrelated four-fold displacement of 0(3) along the < 100 > and Tl along the 
< 110 > directions, and the correlated two-fold displacement of both 0(3) and Tl along the < 110 > directions (with 
the space group Pcca) fll7j , result in essentially the same quality of the fit to the data. In the XAFS data analysis, 
these three models can easily be distinguished as shown above. However, the model we proposed may not be unique. 
More complicated models may generate a similar PPDF to that shown in Fig. 11a. 



B. Comparisons of XAFS data for samples with different T c 

The general picture presented in this paper for the distortions of the TIO layers is similar to that in Ref. 15, but 
with more extensive XAFS data for Tl2Ba2CuOj, and more constraints on possible models. For the XAFS data at 
the Tl edge, both the Tl-Tl and Tl-Ba peaks are split into several peaks which all overlap; at the Ba edge, only the 
Ba-Tl peak is split. Therefore, the XAFS data analysis at the Tl edge is more complicated and technically involved 
than at the Ba edge. Although the Ba edge XAFS data are not quite as sensitive to the Tl distortion as the Tl edge 
XAFS data, they help rule out some possible models which fit the Tl edge XAFS data. 

We also collected Tl Ljjj— edge XAFS data from another sample |fi2ll , Tli.g5Ba2CuO y , with T c =76 K. This sample 
is made in a similar way as described in Sec. IIA. It has a tetragonal structure with lattice constants a — 3.8626 A, 
c = 23.224 A. We also have the Tl L///— edge XAFS data used in the previous study |l5|, which were collected at 
78 K from a tetragonal Tl2Ba2CuOj, sample with T c =89 K and lattice constants a = 3.866 A, c = 23.239 A. All the 
data were reduced in exactly the same way as described in Sec. III. 

The Tl Ljii— edge XAFS data are plotted in Fig. 12 for the three samples with different T c s; solid lines, T c =60 K; 
dotted lines, T c =76 K; dashed lines, T c =89 K. All the data can be fit with the model shown in Fig. 9. The difference 
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in the extracted distances from the constrained fit, among three samples, is within 0.01 A , although the amplitude 
of the Tl-0 peak near 2 A varies significantly among the samples. This clearly shows that the local Tl displacement 
in the TIO layer is not sensitive to the structural symmetry and the transition temperature T c . This local distortion 
is more likely intrinsic to the Tl-2201 compound, which will be further discussed in the following section. 

There is a significant difference in the amplitude of the first peak at r = 2.04 A (c-axis Tl-0(2,3)). The amplitude 
was measured in two different ways. First, the amplitudes were obtained from the constrained fit using the FEFF5 
standards, as described earlier assuming a single peak for the first Tl-0 pair. Second, the data were fit one to another 
in the range of 1.0 ~ 1.8 A, and the relative amplitudes were obtained. The average amplitude is plotted in Fig. 13 
as a function of T c ; it is very close to a linear function. The higher the Tl-0 amplitude of the peak near r = 2 A, the 
lower T c is. This trend is consistent with other experimental results [p H^ffll ; a decrease in oxygen content of about 
0.1^0.15 per formula unit changed T c from to about 80~90 K |4^j4(| although the absolute oxygen content was 
not well determined |||,|5|j4|] . 

The data plotted in Fig. 13 are from three samples with different Tl and O contents which are not precisely 
determined. It has been reported that both the Tl and O contents can change T c significantly pTj. The amplitude 
of the Tl-0 peak plotted in Fig. 13 is from the contributions of the 0(2) and the 0(3) along the c-axis, and possibly 
an 0(4) || if it exists. This amplitude can be affected in several ways: (a) a change in the O content at the 0(2), 
0(3) and 0(4) |] sites, (b) the O movement between the 0(3) and 0(4) sites, (c) changes in the Tl-0(2) and Tl-0(3) 
(along the c-axis) distances, (d) a variation in the Tl content. From the limited data we have, it is not clear what is 
the main reason for the variation of the measured Tl-0 amplitude which seems much too large to be explained solely 
by a variation in O content. However, Fig. 13 clearly shows that T c is very sensitive to the O environment around 
the Tl atom. It is generally believed that the Cu02 layer is the superconducting layer, thus this result also implies 
that a charge transfer occurs between the Cu02 and TIO layers and it is mainly controlled by the O arrangement 
around the Tl, not the Cu substitution for Tl. 



C. The Tl and Cu valences 

The atomic valence (V) can be estimated from the obtained local structure using the Brown- Altermatt-Zachariasen 
formulae which describe the relationship between the bond length (r) and the bond valence (s) 

V = Ssj = H(n/ro)~ N (2) 

or 

V = 'Ss i = Eexp[(r -r i )/B], (3) 

where r a , N and B are empirically determined parameters p8| , ^9| . For Tl 3+ -0, we used Eq. (2) with the most recent 
available parameters, r o =2.04 and N—6A |5(| and the structural model (Fig. 9) constructed based on our XAFS 
results for the sample with the highest O concentration. The calculated Tl valence is (2.95±0.05)+, very close to 
the formal valence, 3+ (Tl2 + Ba2 + Cu 2+ Og _ ). This is also consistent with the photoemission |3l|]52] ] and the x-ray- 
absorption near-edge spectra |H| measurements which indicate a Tl valence slightly less than 3+. If both the Tl and 
0(3) atoms were at their ideal sites (4e), the calculated effective valence for Tl would be 2.77+, which is too small 
compared to the formal valence and not stable for the Tl ion (the Tl ion should be either Tl 1+ or Tl 3+ ). This is 
probably part of the reason why the TIO layer is distorted. 

We also calculated the Cu valence, using Eq. (3) with B = 0.37, r a = 1.679 A jl9| and the local structure 
derived from XAFS (Table I). The calculated Cu valence is (2.15±0.05) + , which is very close to the Cu(2) valence 
in YBa 2 Cu 3 7 (2.20+) ||-|| and Cu in Tl 2 Ba 2 CaCu20 8 (2.11+) || calculated using the same formula. Our 
estimated effective valences for Cu and Tl are qualitatively in agreement with the charge transfer (or local charge) 
picture between the TIO and Cu0 2 layers @,i.e., T1 2 3 + CuO — > [T1 2 3 ]^ + [CuO] 5+ . It is generally believed 
that the extra holes in the Cu0 2 are predominantly at the O site [57|-|59[|. However, it is not clear where the extra 
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electrons are in the TIO layers. We note that the extra holes in the CuC>2 layer are not solely introduced by the 
reduction of the Tl valence. The variation in O content can also change the hole concentration in the CuC>2 layer. 
A decrease in O content of about 0.1 per formula unit, which corresponds to a decrease in hole concentration of 0.2, 
increases T c from ^0 to ~90 K j4^,^(|. However, there is a controversy in the literature about the absolute, optimum 
O content for samples with the highest T c [ 1 1 ll,E5|,E(| . 



D. The question of long range order 

Long range order modulations have been observed in Bi-Sr-Ca-Cu-0 superconductors and related materials p3[ . 
Some long range order in the TIO layers is also observed in Tl-2201 and Tl-2212 by electron diffraction and high- 
resolution transmission electron microscopy ||l8|-f22|. The modulation wave vector is about < —0.16,0.08,1 > for 
tetragonal Tl-2201 (a = b ~ 3.8 A), < 0.08, 0.24, 1 > for orthorhombic Tl-2201 (a ~ b ~5.4 A), and < 0.17, 0, 1 > for 
tetragonal Tl-2212. The superlattice for the first case would have a large unite cell, 6a x 12& in the ab-plane. However, 
these superlattices fl|,|2l[] are usually observed only with selected-area electron diffraction and are much weaker in 
intensity than in the Bi-Sr-Ca-Cu-0 superconductors which exhibit strong supercell diffraction structure |^3|| and 
c-axis displacemnets in the TIO layers are small (see Sec. VI. A). Such superstructures are not always observed in 
the tetragonal Tl-2201 samples [2^ |. Consequently, this long range order of the lattice distortions in the Tl-based 
superconductors seems weak, sample dependent, and may not be intrinsic. Since both the Tl and 0(3) atoms can be 
shifted in two equivalent directions in the model we have proposed, material with this distorted local structure might 
well have a superstructure over a limited range, which is qualitatively consistent with the above electron diffraction 
measurements. There are many ways to construct superstructures with periodicity of 6a x 12b using the distorted 
2a x 2b unit cell shown in Fig. 9 for the TIO layers. Two possible examples that are probably over simplified are 
shown in Fig. 14. Here the Tl displacements (not shown) are the same from one 2a x 2b cell to another, but the 
direction of the O displacements varies. In the upper panel, the O positions are not periodic; some O atoms move 
together or apart, while others move in the same direction. In the lower panel, the direction of the O displacements 
is rotated by 90° in the right half of the supercell. 



VII. CONCLUSIONS 

The local structure of Tl2Ba2CuOj, , with an orthorhombic phase and T c =60 K, has been studied at the Cu K— , 
Ba K— , and Tl L///— edges using the XAFS technique. Several new approaches have been implemented, e.g., direct 
structural simulations up to 7 A using the FEFF5 code, constrained fits, and comparisons of PPDFs for different 
models. 

Both the qualitative and quantitative data analyses clearly show that the Cu, 0(1), 0(2), and Ba atoms are at 
their ideal sites in the unit cell as given by the diffraction studies, while the Tl and 0(3) atoms are displaced from 
the sites suggested by the average crystal structure. The Tl-Tl distance at r=3.5 A between the TIO layers remains, 
but the Tl-Tl distance at 3.9 A in the TIO layer is not observed and the Tl-Ba and Ba-Tl peaks are very broad. The 
shorter Tl-0(3) distance in the TIO layer is about 2.33 A, significantly shorter than the distance calculated with both 
the Tl and 0(3) atoms at their ideal 4e sites ( x = y =0 or j). An excellent fit to the XAFS data can be achieved 
with a correlated displacement model shown in Fig. 9. The fitting results show that the Tl atom is displaced along 
the < 110 > direction from its ideal site by about 0.11 A and the 0(3) atom is shifted from the 4e site by about 0.53 
A roughly along the < 100 > direction. This model also fits very well to the XAFS data collected from two tetragonal 
Tl2Ba2CuO y samples with different T c s. 

The XAFS data do not support the uncorrelated displacement models proposed by diffraction investigations j7j,^],[l0| . 
Our model is similar to that proposed by the PDF analyses of the neutron scattering data on Tl-2212 ||l6fl , but is 
more specific about the relative shift of the Tl atoms between the two consecutive TIO layers. However, the Tl 
displacement, 0.11 A, obtained in our study, is only one third of the value (0.32 A) suggested by the PDF. 
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The estimated Cu and Tl valences imply a charge transfer between the CuC>2 and T10 layers which is consistent 
with other experimental results ]l5] , fi2^|5"l| , [52"| . The estimated Tl valence also supports the distorted structure of the 
T10 layer obtained from our XAFS analysis. The estimated Tl valence (2.95+) in the distorted TIO layer is much 
closer to its formal valence (3+) than that (2.77+) estimated from the ideal TIO layer. 

A comparison of the XAFS data at the Tl Lju— edge from three Tl2Ba2CuOj, samples clearly shows a correlation 
between the transition temperature and the local O environment around the Tl atoms. This implies that the hole 
concentration in the CuC>2 layer is mainly controlled by the O arrangement about the Tl in those samples, and not 
by Cu substitution for the Tl. 
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FIG. 1. The crystal structure for Ti2Ba2CuOi, showing the definitions for the different O sites. Half of the unit cell is shown 
and the whole unit cell is symmetric about the Cu02 plane. The dotted lines indicate the atomic positions in the half unit cell. 
The structural parameters used here are from Ref. 4, in which the 0(3) atom is displaced from the 4e site along the a-axis by 
0.37 A. 

FIG. 2. The k-space XAFS data at the Cu K-, Ba K-, and Tl L///-edges in Tl 2 Ba 2 Cu02 / . 
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FIG. 3. (a-b) The absorptance above the Ba K-edge (a) and the Tl Lj/j-edge (b) in T^I^CuOy; solid lines, experimental 
data; dotted lines, fits to the data using cubic splines; dashed lines, fits to the data using cubic splines plus the extracted back- 
ground from BaO and TI2O3, respectively. The dotted lines and dashed lines are used to simulate the "free-atom absorptance". 
(c-d) The Fourier transform of the Ba K-edge (c) and the Tl L/j_f-edge (d) XAFS data from Ti2Ba2CuO y , using the "free-atom 
absorptance" obtained from cubic spline fits (dotted line) and cubic spline plus corresponding extracted background fits (solid 
lines), respectively. The transform ranges are over 3.5 to 15.5 A -1 at the Ba edge and 3.3 to 15.5 A" 1 at the Tl edge, and both 
are Gaussian broadened by 0.3 A -1 . The fast oscillatory curves are the real parts of the complex transforms and the envelope 
curves are the magnitudes of the transforms. 

FIG. 4. The Fourier transform of the Cu K-edge XAFS, k\(k), to real space for Ti2Ba2CuO y . The transform range is over 
3.1 to 16.5 A -1 , Gaussian broadened by 0.3 A -1 . The fast oscillatory curve is the real part of the complex transform and 
the envelope curve is the magnitude of the transform. The solid lines are the experimental data and the dotted lines are the 
simulations (without a fit) using the FEFF5 code with the average crystallographic structure Q including all the pairs within 
7.5 A. 

FIG. 5. The Fourier transform of the Ba K-edge XAFS, k\(k), to real space for Ti2Ba2CuO y . The transform range is over 
3.5 to 15.5 A -1 , Gaussian broadened by 0.3 A -1 . The fast oscillatory curve is the real part of the complex transform and 
the envelope curve is the magnitude of the transform. The solid lines are the experimental data and the dotted lines are the 
simulations using the FEFF5 code, (a) with the average crystallographic structure [Q , and (b) with the local distorted structural 
model (Fig. 9), including all the pairs within 7.5 A. Part (c) shows a small section of the x— axis enlarged: experimental data 
(solid lines), simulations from part (a) (dotted lines), and simulations from part (b) (dots). 

FIG. 6. The Fourier transform of the Tl Lj/j-edge XAFS, k\(k), to real space for Ti2Ba2CuO y . The transform range is 
over 3.3 to 15.5 A -1 , Gaussian broadened by 0.3 A -1 . The fast oscillatory curve is the real part of the complex transform 
and the envelope curve is the magnitude of the transform. The solid lines are the experimental data and the dotted lines are 
the simulations using the FEFF5 code, (a) with the average crystallographic structure 0], and (b) with the local distorted 
structural model (Fig. 9), including all the pairs within 7.5 A. Note that no fit of the simulated XAFS to the experimental 
data has been made (for more information, see text). 

FIG. 7. The Fourier transform of k\(k) at the Cu K-edge for Tl2Ba2CuO y . The transform range is the same as in Fig. 4. 
Both the real part and the magnitude of the transform are plotted. The solid lines are the experimental data and the dotted 
lines are the fits to the data. The fitting range is 1.0 to 3.9 A. 

FIG. 8. The Fourier transform of k\(k) at the Ba K-edge for Ti2Ba2CuO y . The transform range is the same as in Fig. 5. 
Both the real part and the magnitude of the transform are plotted, (a) The comparison between the experimental curves and 
fits to the data. The fitting range is 1.2 to 3.8 A. The solid lines are the experimental data and the dotted lines are the fits, 
(b) The contributions of the Ba-O, -Cu, and -Ba peaks in the fit. (c) The contributions of the Ba-Tl pairs. The top trace is 
the sum of the Ba-Tl peaks, followed by the three individual peaks. 

FIG. 9. The projection of the local structural model in the ab-plane for the distorted TIO layers. Two successive layers are 
shown. The best overall agreement between the model and the XAFS results is achieved when the Tl atom is shifted along the 
< 110 > direction by 0.11 A and the O atom in the TIO layer is shifted along the [100] and/or [010] directions by about 0.53 
A from its ideal site. The O displacement along the [010] direction is not shown in the figure. Note that the displacements of 
the Tl and O atoms are correlated, not random. 
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FIG. 10. The Fourier transform of k\{k) at the Tl L///-edge for Tl 2 Ba2CuO y . The transform range is the same as in Fig. 
6. Both the real part and the magnitude of the transform are plotted, (a) The comparison between the experiment results and 
fits to the data. The fitting range is 1.2 to 4.2 A. The solid lines are the experimental data and the dotted lines are the fits 
to the data, (b) The contributions of the Tl-0 peaks in the fit. (c) The contributions of the Tl-Tl pairs. The top trace is the 
sum of the Tl-Tl peaks, followed by the three individual peaks, (d) The contributions of the Tl-Ba pairs. The top trace is the 
sum of the Tl-Ba peaks, followed by the three individual peaks. 

FIG. 11. The Tl-Tl pair-distribution-function (PDF) from (a) the correlated Tl displacement model (refer to Fig. 9), (b) the 
ideal structure, (c) the uncorrelated Tl displacements along the < 110 > direction, and (d) the uncorrelated Tl displacements 
along the < 100 > direction. The dotted lines are the starting PDFs and solids are the extracted Tl-Tl PDFs from XAFS using 
the corresponding structural models and the constrained fits. For more information, see Sec. VI. A. 

FIG. 12. The Fourier transform of the Tl L///-edge XAFS, k\{k), to real space. The transform range is over 3.2 to 15.5 
A -1 , Gaussian broadened by 0.3 A -1 . Both the real part and the magnitude of the Fourier transform are plotted. The data 
are collected at T=80±2 K from three samples with T c =60 K (orthorhombic, solid lines), T c =76 K (tetragonal, dotted lines), 
and T c =89 K (tetragonal, dashed lines), respectively. For more information, refer to Sec. VLB. 

FIG. 13. The amplitude of the nearest Tl-O peak near r — 2 A versus T c for three samples presented in Fig. 12. 

FIG. 14. Two possible models for the 6a x 126 superlattice. The Tl displacements are the same from one 2a x 2b unit cell 
to another while the O displacements vary. 
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TABLE I. Cu local structure obtained from XAFS and diffraction, "nbrs" indicates the weighted number of neighbors. The 
estimated errors are: number of neighbors, ± 15 %; distances, ± 0.02 A. 



Pairs 


r(A) 


XAFS 
nbrs 


(7(A) 


Diffraction 

r(A) 


nbrs 


Cu-O(l) 


1.93 
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1.93 
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Cu-0(2) 
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0.093 


2.72 
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Cu-Ba 


3.35 
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0.053 


3.35 


8 


Cu-Cu 


3.86 


3.2 


0.046 


3.87 


4 


Cu-O(l) 


4.31 


7.7 


0.056 


4.32 


8 



TABLE II. Ba local structure from XAFS, diffraction, and the model. The values listed in the brackets are from other 
diffraction studies [p|-[il"[. "nbrs" indicates the weighted number of neighbors. In the model, the Tl atoms are shifted from the 
4e site along the < 110 > direction by 0.11 A and the 0(3) atoms are shifted from the 4e site by 0.53 A along the < 100 > 
direction. The estimated errors are: number of neig hbors, ± 15 %; distances, ± 0.02 A. 



Pairs 




XAFS 




Diffraction 




Model 






r(A) 


nbrs 


a(A) 


r(A) 


nbrs 


r(A) 


nbrs 


Ba-O(l) 


2.73 


4.0 
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TABLE III. 


Tl local structure obtained from XAFS, diffraction, and the model. 
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the < 110 > direction by 0.11 A and the 0(3) atoms are shifted from the 4e site by 0.53 A alon; 


g the 


< 100 > direction. The estimated errors are: number of neighbors, ± 15 %; distances, ± 0.02 A. 
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